arXiv:1503.01860vl [astro-ph.CO] 6 Mar 2015 


Draft version March 9, 2015 

Preprint typeset using ET^X style emulateapj v. 5/2/11 


USING MEMBER GALAXY LUMINOSITIES AS HALO MASS PROXIES OF GALAXY GROUPS 

Yi Lu\ Xiaohu Yang 1 2 3 ' 2,3 , Shiyin Shen 1 
Draft version March 9, 2015 

ABSTRACT 

Reliable halo mass estimation for a given galaxy system plays an important role both in cosmology and 
galaxy formation studies. Here we set out to find the way that can improve the halo mass estimation for those 
galaxy systems with limited brightest member galaxies been observed. Using four mock galaxy samples con¬ 
structed from semi-analytical formation models, the subhalo abundance matching method and the conditional 
luminosity functions, respectively, we find that the luminosity gap between the brightest and the subsequent 
brightest member galaxies in a halo (group) can be used to significantly reduce the scatter in the halo mass 
estimation based on the luminosity of the brightest galaxy alone. Tests show that these corrections can signif¬ 
icantly reduce the scatter in the halo mass estimations by ~ 50% to ~ 70% in massive halos depending on 
which member galaxies are considered. Comparing to the traditional ranking method, we find that this method 
works better for groups with less than five members, or in observations with very bright magnitude cut. 
Subject headings: large-scale structure of universe - dark matter - galaxies: halos - methods: statistical 


1. INTRODUCTION 

In the current scenario of galaxy formation, galaxies are 
thought to be formed and reside in cold dark matter haloes. 
Studying the galaxy-halo connection observationally provides 
one with important information about the underlying pro¬ 
cesses in galaxy formation and evolution. In recent years 
great progress has been made in establishing the halo-galaxy 
link via the so called halo occupation models (Jing, Mo & 
Borner 1998; Mo, Mao & White 1999; Peacock & Smith 
2000; Zheng et al.2007) and the conditional luminosity func¬ 
tions (Yang et al. 2003; van den Bosch et al. 2003; Cooray et 
al. 2006). This halo-galaxy connection can also be made di¬ 
rectly using galaxy groups which are defined as sets of galax¬ 
ies that reside in the same dark matter halo, as is done in Yang 
et al. (2005c,d), Collister & Lahav (2005), van den Bosch et 
al. (2005), Robotham (2006), Zandivarez et al. (2006), Wein- 
mann et al. (2006a,b). 

Numerous group catalogues have been constructed from 
galaxy redshift surveys, including the 2-degree Field Galaxy 
Redshift Survey (2dFGRS) (Eke et al. 2004; Yang et 
al. 2005a), the high-redshift DEEP2 survey (e.g.. Crook et 
al.2007); and the Sloan Digital Sky Survey (SDSS) (e.g., 
Berlind et al. 2006; Yang et al. 2007; Tago et al. 2010; Nurmi 
et al. 2013), the Galaxy and Mass Assembly (GAMA) ob¬ 
servations (e.g. Robotham et al. 2011) using different meth¬ 
ods, ranging from the traditional friends-of-friends (FOF) al¬ 
gorithm, to the hybrid matched filter method (Kim 2002) and 
the “MaxBCG” method (Koester et al. 2007). Relevant to 
the present paper is the halo-based group finder developed by 
Yang et al. (2005a), which groups galaxies according to their 
common halos expected from the current CDM model. This 
group finder is suited to study the relation between galaxies 
and dark matter haloes over a wide dynamic range in halo 
mass, from rich clusters of galaxies to poor groups of galax¬ 
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ies, as tested with mock galaxy surveys, and has been applied 
to 2dFGRS (Yang et al. 2005a), SDSS DR2 (Weinmann et al. 
2006a), DR4 and DR7 (Yang et al. 2007). 

One of the key steps in the halo-based group finder is the 
estimation of halo masses of candidate galaxy groups. In gen¬ 
eral, the halo masses of a group can be estimated based on the 
velocity dispersion of their member galaxies. However, ex¬ 
cept for a small fraction of very rich groups with a large num¬ 
ber of satellite galaxies (e.g., Carlberg et al. 1996; 1997), this 
method is not very reliable for small groups with only a few 
members. In general, stacking of satellite galaxies for given 
central galaxies are used to obtain their halo masses (e.g., Er¬ 
ickson et al. 1987; Zaritsky et al. 1993; McKay et al. 2002; 
Brainerd & Specian 2003; Prada et al. 2003; van den Bosch 
et al. 2004; Becker et al. 2007; Conroy et al. 2007; Norberg 
et al. 2008; More et al. 2009b, 2011; Li et al. 2012). 

Group total luminosity (e.g. Yang et al. 2005a; 2007) or 
richness (e.g. Berlind et al. 2006; Becker et al. 2007; An- 
dreon & Hurn 2010; Hao et al. 2010) can also be used as 
halo mass indicator. As demonstrated in Yang et al. (2007), 
the halo mass is tightly correlated with the total luminosity of 
member galaxies. Unfortunately not all member galaxies can 
be observed because galaxy samples are usually magnitude 
limited. In practice, one can only measure a characteristic 
group luminosity, L G , which is defined as the total luminosity 
of member galaxies brighter than a given luminosity thresh¬ 
old. In the SDSS group catalogue of Yang et al., this luminos¬ 
ity threshold is chosen to be 01 M r — 5 log h = —19.5, so that 
one can reach to a depth of z = 0.09. Here 01 M T . — 5 log h 
is the absolute r-magnitude, K- and E-corrected to z = 0.1, 
the typical redshift of galaxies in the SDSS redshift sample. 
Assuming that L G increases monotonically with halo mass 
Mh, one can then obtain the halo masses of galaxy groups by 
abundance matching between groups and dark matter halos, 

/•oo r-oo 

/ n G (L' G )dL' G = n h (M' h )dM ' h , (1) 

J Lg J Mh 

where n G {L G ) is the number density of groups at character¬ 
istic luminosity L G , and rih{Mh) is the halo mass function 
for the cosmology adopted. In Yang et al. (2007), Mh — L G 
relation so obtained is used iteratively in the group finder to 
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link galaxies in the same dark matter halo, and to assign halo 
masses to individual galaxy groups. 

As shown in Yang et al. (2005a; 2007) using mock galaxy 
catalogs, the halo based group finder and the mass assignment 
scheme work well for a moderately deep and large survey, 
such as the 2dFGRS and SDSS. However, the method may 
not work as well for shallow surveys, such as the 6dFGRS 
(Jones et al. 2009), and for groups selected from high redshift 
surveys, such as DEEP2 (Willmer et al. 2006), COSMOS (II- 
bert et al. 2009) and GAMA (Driver et al. 2011) where only 
a few (in most cases one or two) brightest member galaxies 
are observed in each dark matter halo. In case the survey 
volume is very small or difficult to calculate because of the 
bad survey geometry, the halo mass estimation based on the 
ranking of Lq may become unachievable. Moreover, there 
are also some catalogues which only consists of a certain type 
of galaxies, for example, the SDSS-III’s Baryon Oscillation 
Spectroscopic Survey (BOSS) (Dawson et al. 2013) which 
map the spatial distribution of luminous red galaxies (LRGs). 
It is unclear if the method can work well for those catalogues. 
Thus, an alternative estimate of the halo masses is required 
for such groups. 

To estimate the halo masses for galaxy systems with only 
one or two brightest member galaxies, one may make use of 
the central-host halo relation obtained from conditional lumi¬ 
nosity function (CLF; e.g. Yang et al. 2003; van den Bosch 
et al. 2003), the subhalo abundance matching (SHAM; e.g.. 
Vale & Ostriker 2006), or from predictions of semi-analytical 
model (SAM; e.g., Kang et al. 2005). As shown in Yang 
et al. (2008) and Cacciato et al. (2009), the typical scatter 
in logL c for halos of a given mass is about 0.15. However 
for massive halos L c oc M ^ 0 25 , and so the scatter in halo 
mass for a given L c can be substantial. Thus the central (or 
the brightest) galaxy alone can not provide a reliable estima¬ 
tion of the halo mass, especially for massive halos. Recently, 
suggestions have been made to the magnitude gap between 
the brightest and second brightest galaxies as another param¬ 
eter to describe the halo mass in addition to L c (e.g., Hearin 
et al. 2013a; More 2012; Shen et al. 2014). In this pa¬ 
per, we investigate how halo masses are correlated with group 
properties, such as the central galaxjQ luminosity, the satellite 
member galaxy luminosity or the luminosity gap, defined as 
log Tgap = log L c — log L s , where L s is the luminosity of the 
i-th brightest member galaxies. For example, if a group only 
has two members, then L s is the luminosity of the second 
brightest galaxy. If a group contains five member galaxies, 
we choose to use the luminosity of the fifth brightest member 
galaxy. The goal of this paper is to improve the halo mass 
estimation using both L c and the luminosity gap. The investi¬ 
gations are carried out based on four mock samples generated 
using SAM, SHAM, and two CLF models, respectively. 

The outline of this paper is as follows, ^describes the four 
mock samples used in this paper that are constructed using the 
SAM, SHAM, and CLF, respectively. In Sj3]we present the re¬ 
lationships between the hal o ma ss and the total group (or cen¬ 
tral galaxy) luminosity. In T3.3I we show how the luminosity 
gap affects the L c — Mh relation and can be used to reduce 
the scatter within this relation. In 113.41 we investigate the im¬ 
provement of this reduction by involving more fainter mem¬ 
ber galaxies in calculating luminosity gap. In 10] we compare 
the performance of two halo mass estimation methods, which 

4 Throughout this paper, we refer to the brightest galaxy in a halo as the 
central galaxy. 


are based on luminosity gap (refer to ‘GAP’) and traditional 
ranking method (‘RANK’) under different circumstances, re¬ 
spectively. Finally, we summarize our results in 1|5] 

2. SAMPLES OF MOCK GROUPS 

In this paper, we make use of four sets of mock galaxy cat¬ 
alogs. The first is constructed by Guo et al. (2011; here¬ 
after Gil) using a semi-analytic model of galaxy formation 
applied to dark matter halo merging trees obtained from the 
‘Millennium Run’ N-body simulation (Springel et al. 2005). 
The cosmological parameters adopted in the simulation are 
Q m = 0.25, Oa = 0.75, h = 0.73 and a CDM spectrum 
with an amplitude specified by erg = 0.9. This set of parame¬ 
ters is different from that obtained in recent observations (e.g. 
Planck Collaboration et al. 2013; Mandelbaum 2013). How¬ 
ever, for our test of accuracy of halo mass assignment, this 
is not a big issue. The simulation was performed with GAD- 
GET2 (Springel 2005) using 2160 3 dark matter particles in a 
periodic cubic box with a side length Lbox = 500/i _1 Mpc (in 
comoving units). The mass of a particle is 8.6 x 10 8 h. _1 MQ. 
In their modeling, Gil included the growth and activity of 
the central black holes, as well as their effects on suppress¬ 
ing the cooling and star formation in dark matter halos. Since 
we are interested in finding an accurate halo-mass proxy from 
observed luminosities, here we only use the halo mass and r- 
band luminosity of each galaxy given by Gil. We refer the 
reader to Gil for the details of the semi-analytic model (see 
also in Guo et al. 2013). In this paper, we only make use of 
these galaxies with luminosity log L > 8.0. This set of mock 
galaxy catalog is referred to as ‘SAM’. 

The second mock galaxy catalog is constructed by Hearin 
et al. (2013b). They used the Bolshoi N-body simulation 
(Klypin et al. 2011) which models the cosmological growth 
of structure in a cubic volume with side length Lbox = 
250 Ji -1 Mpc within a standard ACDM cosmology with pa¬ 
rameters: f l m = 0.27, Ha = 0.73, h = 0.7, n = 0.95 
and (xg = 0.82. After utilizing the ROCKSTAR (Behroozi et 
al. 2012) halo finder to identify halos and subhalos, subhalo 
abundance matching (SHAM) technique was adopt to asso¬ 
ciate galaxies with dark matter halos. This technique, called 
as age distribution matching, is a two-phase algorithm for as¬ 
signing both luminosity and color to the galaxies located at 
the center of halos (see Hearin et al. 2013a; 2013b for de¬ 
tails). Here, we make use of their galaxy catalog in the cubic 
box that contains 244784 galaxies at redshift 2 = 0.1 with 
r-band absolute magnitude 01 M r — 51og/i < —19.0. 

The third mock galaxy sample used in this paper is con¬ 
structed by populating dark matter halos obtained from nu¬ 
merical simulations with galaxies of different luminosities ac¬ 
cording to the conditional luminosity function (CLF; Yang et 
al. 2003). The algorithm of populating galaxies is similar to 
that outlined in Yang et al. (2004) but here updated to the 
CLF in r-band. The CLF, f l»(L | A7/,), is defined to be the aver¬ 
age number of galaxies of luminosities L ± dL/2 that reside 
in a halo of mass Mh (see Yang et al. 2003). Here we give 
a brief description of the algorithm of assigning luminosity 
to each member galaxy. First, we write the total CLF as the 
sum of the CLFs of central and satellite galaxies (Yang et al. 
2003; van den Bosch et al. 2003; Vale & Ostriker 2004, 2006; 
Cooray 2006; van den Bosch et al. 2007; Yang et al. 2008): 

*(L\M h ) = 3> cen (L\M h ) + $ sat (L|A4). (2) 

We assume the contribution from the central galaxies to be a 
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lognormal distribution: 


^cen {L\M h ) dlog L 
1 

= —j =—exp 
v 27 TCT c 


(logi - log L c ) 2 


(3) 


d log L , 


where a c is a free parameter, which expresses the scatter in 
log L of central galaxies at fixed halo mass, and log L c is the 
expectation value for the ( 10 -based) logarithm of the luminos¬ 
ity of the central galaxy. For the contribution from the satellite 
galaxies we adopt a modified Schechter function which de¬ 
creases faster at the bright end than the Schechter function: 


$ sat (L|M) dlog L 



(«.+i) 

exp 


(4) 



ln( 10 ) dlogL. 


The parameters L c , a c , (/>*, a s and L* are all functions of the 
halo mass Mh- 

Following Yang et al. (2008) and Cacciato et al. (2009), we 
assume that a c is a constant independent of halo mass, and 
that the L c — Mh relation has the following functional form. 


L c (Mh) — L 0 


(1 + Affc/MiJ-ri-T* ' 


(5) 


This model contains four free parameters: a normalized lumi¬ 
nosity, Lq, a characteristic halo mass, Mi, and two slopes, 71 
and 72. For satellite galaxies we use 


log L*(M h ) = log L c (M h ) - 0.25 , ( 6 ) 

a s {M h )=a s (7) 

(i.e., the faint-end slope of $ sat (L| A7/, ) is independent of halo 
mass), and 

log \4>* s {M h )\ = b 0 + 61 (log M 12 ) + 6 2 (logMi 2 ) 2 , ( 8 ) 


the mean obtained in step ( 1 ) to obtain the number of satel¬ 
lite galaxies; (3) assign a luminosity to each of these satellite 
galaxy according to Eq. Q. Note that satellite galaxies are 
allowed to be brighter than the central galaxy. However, we 
note that in our subsequent analysis, following the observa¬ 
tional definition, the central galaxy is defined as the brightest 
galaxy in a halo. We refer to this set of mock galaxy catalog 
as ‘CLF1’. 

Finally, in order to test if our analysis is cosmological de¬ 
pendent or not, we construct another set of mock galaxy cat¬ 
alog which we call as ‘CLF2’. The cosmological parameters 
adopted for this set of mock galaxy catalog are those given by 
WMAP3: Q m = 0.238, fl A = 0.762, h = 0.73, n = 0.95 
and erg = 0.75. We adopt the set of the best fit CLF param¬ 
eters that are consistent with WMAP3 cosmology, provided 
in Cacciato et al. (2009), with (logMi = 11.07, logLo = 
9.935 ,71 = 3.273 ,72 = 0.255, a c = 0.143) (see Cacciato et 
al. 2009 for more details). Here again, we use the halo cat¬ 
alog from the ‘Millennium Run’ N-body simulation to con¬ 
struct the mock catalog. Since the cosmology adopted in the 
Millennium simulation is different from the WMAP3 cosmol¬ 
ogy, we scale the halo masses in the Millennium simulation to 
the WMAP3 cosmology using abundance matching according 
to the mass functions of WMAP3 cosmology and the Millen¬ 
nium simulation. Such abundance matching ignores the dif¬ 
ference in the spatial correlations between halos in different 
cosmologies, but should be valid for our analysis which is 
based only on galaxy occupation in individual halos. 

Among the above four mock catalogs, although the SAM 
has more ‘physics’ in it, but doesn’t really fit the data (e.g., 
see Weinmann et al. 2006; Liu et al. 2010; Nurmi et al. 
2013), whereas the CLF1, CLF2 and SHAM mocks are not 
based on physical models of galaxy formation but do fit the 
observational data. Here in our subsequent probes, we make 
use of CLF1 as our fiducial sample, and take others mainly 
for comparisons. 


with Mi 2 = Mh/{ 10 12 6. -1 Mq). This set of CLF parameter¬ 
ization thus has a total of nine free parameters, characterized 
by the vector 

A CLF = (logMi,logL 0 ,7 i, 72,CT c ,a s ,60,61,62) • ( 9 ) 

The above functional forms do not have ample physical mo¬ 
tivations but were found to adequately describe the observa¬ 
tional results obtained by Yang et al. (2008) from the SDSS 
galaxy group catalog. They were adopted in Cacciato et al. 
(2009) to make the model prediction for the galaxy-galaxy 
lensing signals, and very recently in van den Bosch et al. 
(2013), More et al. (2013) and Cacciato et al. (2013) to 
constrain both the CLF and cosmological parameters with the 
SDSS clustering and weak lensing measurements. Here we 
adopt the set of the best fit CLF parameters that are obtained 
for the Millennium simulation cosmology, with (log Mi = 
10.954, logLo = 9.896,7 i = 5.192, 72 = 0.2415, a c = 
0.1501, a s = -0.6828,6 0 = -0.1611 ,61 = 0.7945, 6 2 = 
-0.1573). 

Galaxies with luminosities \og(L / h~ 2 L Q ) > 8.0 are pop¬ 
ulated in the dark matter halos extracted from the Millen¬ 
nium Simulation. In practice, each halo is assigned a central 
galaxy with a median luminosity specified by Eq. ([5} and a 
log-normal dispersion a c = 0.1501. The satellite galaxies are 
populated via following steps: ( 1 ) obtain the mean number of 
satellite galaxies according to the integration of Eq.® with 
luminosities logL > 8.0; (2) use a Possion distribution with 


3. HALO MASSES FROM LUMINOSITIES OF MEMBER 
GALAXIES 

In order to carry out cosmological and structure forma¬ 
tion investigations using observations of galaxy systems, one 
of the key steps is to obtain reliable estimations of the halo 
masses of the systems selected. The kinematics of member 
galaxies is one of the most commonly used methods to esti¬ 
mate masses of optically selected groups/clusters (e.g., van 
den Bosch et al. 2004; Kargert et al. 2004; More et al. 
2011). However, this method may be significantly affected 
by anisotropies of the velocity dispersion (e.g. Biviano et al, 
2006; White et al. 2010). Furthermore, for groups containing 
only a few members, or when only a few group members can 
be observed, the estimation of velocity dispersion becomes 
noisy, making the dynamical mass estimation unreliable. 

Recently, another halo mass indicator has been used to es¬ 
timate halo masses for galaxy systems over a large dynami¬ 
cal range. This method uses the richness or the characteris¬ 
tic group luminosity as the proxy of halo mass, and so can 
readily be applied once galaxy systems are selected from a 
redshift survey. Richness is one of the easiest quantities of a 
galaxy system to observe, but there are different ways to de¬ 
fine the richness. For instance, one can define it as the number 
of member galaxies within the red sequence or above a cer¬ 
tain luminosity threshold. Gladders & Yee (2000, 2005) used 
the red sequence richness to measure the mass of the halos 
from photometric data. The characteristic group luminosity 
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log [L 0 / h _ 2 L s ] log (L 0 / h - 2 L e ] 

Fig. 1.— The halo mass M v.s. characteristic group luminosity Lq in the CLF1 (top left panel), CLF2 (top right panel), SAM (bottom left panel) and SHAM 
(bottom right panel) mock catalogs. Here results are shown for the median (symbols) and 68% confidence level (error bars) of log Mh within each log Lq bin. 
The vertical green dashed line in each panel shows the luminosity limit at 01 M r — 5 log h = —19.5. 


(or stellar mass), which is defined as the total group luminos¬ 
ity (stellar mass) of member galaxies above a certain lumi¬ 
nosity threshold, is presumably a better mass indicator than 
the richness. This is especially the case for very poor galaxy 
systems dominated by a single luminous galaxy, where the 
galaxy luminosity is known to change as a function of halo 
mass over a large range. 

3.1. Using the characteristic group luminosity 

As demonstrated in Yang et al. (2005a; 2007) using mock 
redshift surveys, the total luminosity of all member galaxies 
in a group is tightly correlated with halo mass. This method 
needs the determination of the total group luminosity. How¬ 
ever, a limitation of this method is the need of corrections 
in real observations where only bright galaxies are observed. 
Although one may make corrections using a luminosity func¬ 
tion (e.g. that for the total galaxy population), detailed CLF 
modeling shows that the correction depends significantly on 
halo mass (Yang et al. 2003; van den Bosch et al. 2007; 
Yang et al. 2008). To avoid this uncertainty and as a com¬ 
promise with observational magnitude limit, one may use a 
characteristic group luminosity, defined as the total luminos¬ 
ity of member galaxies above a luminosity threshold. For in¬ 
stance, Yang et al. (2005a; 2007) used a characteristic group 
luminosity, Lq, which is defined as the total luminosity of all 
group members brighter than — 5 log h = —18 for the 

2dFGRS group catalogue, or of all group members brighter 
than 01 M r — 5 log ft, = —19.5 for the SDSS DR7 group cat¬ 
alogue, as the halo mass proxy. Tests using mock galaxy and 
group catalogues show that the characteristic group luminos¬ 
ity is tightly correlated with the mass of the dark matter halo 
hosting the group (Yang et al. 2005a). Once the character¬ 


istic group luminosities of all the groups are obtained, one 
can assign groups with halo masses according to abundance 
matching (e.g. Mo, Mao & White 1999) between the group 
luminosity function and the halo mass function predicted with 
a given cosmology (see Yang et al. 2007 for details). 

Fig. [I] shows the relation between the characteristic group 
luminosity Lq (the total luminosity of all group members 
brighter than 01 M r . —5 log ft = —19.5) and the halo mass Mh 
obtained directly from the CLF1, CLF2, SAM and SHAM 
mock galaxy catalogs in different panels as indicated, respec¬ 
tively. The symbols in each panel are the median values while 
the error-bars are the l-er variation ( 68 % confidence level). 
There is a clear tight correlation between the characteristic 
luminosity and the halo mass. The typical l-er variations in 
log Mh given by the SAM sample are ~ 0.2 for massive ha¬ 
los and ~ 0.4 for intermediate mass {Mh ~ 10 12 ' 5 ft _ 1 M 0 ) 
halos. The corresponding variations given by the CLF1, 
CLF2 and SHAM samples are ~ 0.1 for massive halos and 
~ 0.3 in the intermediate mass range. Note that, the very 
asymmetric error bars at low luminosity end for SAM sam¬ 
ple are caused by a combination of the color modeling of 
central galaxies in halos around 10 118 ft _1 M 0 (see Fig. 2 
for such a feature for central galaxies), the luminosity cut at 
°- 1 M r — 5 log ft = —19.5 and a smaller bin size used here. 

3.2. Using the brightest (central) galaxy 

In the local Universe, the characteristic group luminosi¬ 
ties can be determined for groups down to a halo mass of 
~ 10 12 ft _ 1 M 0 with surveys like the 2dFGRS and the SDSS. 
However, if one goes to higher redshift or uses much shal¬ 
lower observations (e.g. 6 dFGRS), where only a few brightest 
members a group can be observed, corrections for the missing 
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Fig. 2.— The median halo mass, Mh, as function of central galaxy luminosity, L c . The red solid circles show the median values of Mh in given central galaxy 
luminosity bins. Error bars indicate the 68% confidence level around the median values. Here results are shown for mock samples CLF1 (top left), CLF2 (top 
right), SAM (middle left) and SHAM (middle right), respectively. The solid line in each panel is the best fit Mh — L c relation. As a comparison, we also show 
using blue open circles the median log L c in different log M ^ bins with horizontal error bars indicating the 68% confidence levels around the median. The dotted 
line in each panel is the best fit L c — M ^ relation. The bottom two panels show all the fitting lines of L c — Mh relation (bottom left) and M ^ — L c relation 
(bottom right) for four mock samples respectively. 


members have to be made to obtain the characteristic group 
luminosity, making the method unreliable. On the other hand, 
investigations with the CLF model and observed group and 
cluster catalogues have shown that the luminosity of the cen¬ 
tral galaxy L c of a group is correlated with its halo mass Mh 
(e.g., Yang et al. 2003; van den Bosch et al. 2003; Yang et al. 
2007). As an illustration, FigU shows the model predictions 
of the L c — Mh relation from the four mock samples described 
in Section [2] as marked in each panel, respectively. In each 
panel, the blue open circles are the median values of the cen¬ 
tral galaxy luminosity in different halo mass bins, and the hor¬ 
izontal error bars show the 68 % confidence level around the 
median luminosities. As one can see, the median luminosity 
of central galaxy increases monotonically with halo mass, and 
that the amounts of scatter are comparable among four mock 
catalogs. For reference, we fit the L c — Mh relations using the 
functional form given by Eq. (0. The fitting parameters are 
listed in Table [T] and the best fitting results are shown in Fig. 
[2] with the dotted lines. For the CLF1 and CLF2 samples, the 
parameters are adopted the same as those listed in Section [2] 
These fitting formulae describe the median L c — Mh relations 


TABLE 1 

Parameters for the L c — M h relation 


Sample 

log L 0 

log Mi 

71 

72 

CLFl 

9.896 

10.954 

5.192 

0.2415 

CLF2 

9.935 

11.07 

3.273 

0.255 

SAM 

10.02±0.25 

10.95i0.48 

2.59i0.40 

0.36i0.03 

SHAM 

10.25i0.16 

12.94i0.56 

0.36i0.03 

0.22i0.03 


remarkably well. For the most massive halos, the SAM pre¬ 
dicts central galaxies are significantly brighter than the CLF 
model. However, as compared among SAMs and SDSS ob¬ 
servations in G11, the stellar mass functions at massive end of 
SAMs do not show very significant deviation, the above dif¬ 
ference may due to the color modeling of these galaxies. Their 
model central galaxies, especially in ~ 10 118 /t _1 M 0 halos, 
tend to be too blue and bright. This is the case of the wiggle in 
their Mh — L c relation and the resulting much smaller median 
L c — Mh relation. 

Although the L c — Mh relations presented above are useful 
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Fig. 3.— Similar to the results shown in Fig. [2] but here for groups within three different luminosity gap ranges: 0.5 > logLgap > 0.0 (solid circles), 
1.0 > log Lgap > 0.5 (open circles) and logL gap > 1.0 (solid squares). Here again, symbols and error bars represent the median and 68% confidence levels 
of log M in different log L c bins. 


in predicting the luminosities of central galaxies in halos of 
given masses, these relations are not appropriate for estimat¬ 
ing halo mass from L c . Because the luminosity function is 
steep at the bright end, and because the scatter in the L c {Mh) 
relations is quite large, M), will be over-estimated from the 
L c (Mh) relation due to a Malmquist-like bias. In order to 
have an unbiased result, we need to obtain the Mh(L c ) rela¬ 
tion directly instead of using the inverse of the L c (Mh ) rela¬ 
tion. In Fig. [2] we plot the median (red solid circles) and the 
68 % confidence levels (error bars) of halo masses, log Mh, 
as a function of central galaxy luminosity log L c , obtained di¬ 
rectly from the CLF1, CLF2, SAM and SHAM samples. Here 
again the median halo mass increases monotonically with the 
central galaxy luminosity. For low-mass halos, the median of 
the Mh(L c ) relation is similar to that of the L c (Mh) relation, 
but for halos with mass Mh > 10 12 /i _1 M 0 these two rela¬ 
tions are quite different, exactly because of the Malmquist- 
like bias. The scatter in log Mh for a given luminosity in¬ 
creases with log L c , and at the bright end is much larger than 
that in the L c (Mh) relation, particularly in the SAM predic¬ 
tion, clearly due to the shallow slope in the L c (Mh) relation 
(for example, L c oc M^ 0 ' 25 at the massive end for the CLF 
sample) (see More et al. 2009a; 2009b for more detailed dis¬ 
cussions). On the other hand, as pointed out in More et al. 
(2009b), the huge scatter in Mh at large L c seen in SAM is 
inconsistent with the constraints obtained using satellite kine¬ 
matics. We fit the median Mh(L c ) relations with the follow¬ 
ing functional form: 

log M h = exp (log L c - log M a ) ± log M b . (10) 

The best fitting parameters of [log M a , log M b \ are: [9.61 ± 


0.01,10.60±0.01] for the CLF1 sample, [9.69±0.02,10.56± 
0.03] for the CLF2 sample, [10.32 ± 0.03,10.93 ± 0.04] for 
SAM sample and [9.68 ±0.02,10.64 ± 0.03] for SHAM sam¬ 
ple, respectively. The resulting best fit median relations are 
shown in Fig. 2 as the solid lines. Note that the above fitting 
formula are obtained within the halo mass range shown in the 
plot, extrapolation to much larger halo mass is not valid. For 
clarity, we put all the fitting lines of four mock samples in 
the bottom two panels for L c — Mh relations on the right and 
Mh — L c relations on the left. Obviously, for SAM, both rela¬ 
tions are inconsistent with the other three mock samples. As 
we have mentioned before, it’s because the galaxy luminos¬ 
ity in SAM model is overestimated compare to the other three 
mock samples. 

The above results show that the scatter in the Mh(L c ) is 
quite large at the massive end. It is therefore not a good choice 
to use L c alone as a halo mass proxy. On the other hand, we 
see in Fig. 1 that the scatter shown for a given characteristic 
group luminosity is much smaller, which indicates that us¬ 
ing additional member galaxies can give better constrain of 
the halo masses. In what follows we investigate how to get 
a better halo mass proxy by using the luminosities of other 
brightest member galaxies. 

3.3. Using the brightest two galaxies 

To start with, we introduce a ‘luminosity gap’, defined as 
the luminosity ratio between the central and a satellite galaxy 
in the same dark matter halo, L gap = L c /L s . In this section, 
we focus on the brightest satellite galaxy, so that L s = ± 2 , 
where ±2 is the luminosity of the brightest satellite or the sec¬ 
ond brightest among all the member galaxies. Recent studies 
have shown that this luminosity gap may contain important 
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Fig. 4.— The halo mass difference, A log Mh , between the best lit Mh — L c relation and data as a function of luminosity gap log L ga p. The results shown in 
different panels are for central galaxies within different luminosity bins: log L c = 9.9 ± 0.05 (upper-left), 10.2 ± 0.05 (upper-right), 10.5 ± 0.05 (lower-left) 
and 10.8 dz 0.05 (lower-right), respectively. In each panel, red solid circles are results obtained from the CLF1 sample, where the symbols and error bars represent 
the median and 68% confidence levels of A log Black solid line is the best fitting results to the data from CLF1 sample. The results obtained from CLF2, 
SAM and SHAM samples are also shown in each panel for comparison using different symbols as indicated. 


information about halo mass (e.g., Hearin et al. 2013; More 
2012). As found in Shen et al. (2014), L 2 is also correlated 
with group richness which, in turn, is correlated with halo 
mass. Thus a combination of the luminosities of the bright¬ 
est and second brightest galaxies in a halo is likely to provide 
more information regarding the halo mass than the brightest 
galaxy alone. 

The luminosity gap L gap can be measured straightfor¬ 
wardly for groups/halos that have at least two members. 
Within our four samples, the distribution of the luminosity 
gap, log Lg a p, spreads roughly in the range 0.0-3.0. We 
split each of our mock group catalogues into three subsam¬ 
ples according to the value of logL gap : small gap groups 
with 0.5 > logLgap > 0.0; median gap groups with 1.0 > 
log Lgap > 0.5; and large gap groups with logL gap > 1.0. 
Fig. [3] shows the same Mh(L c ) relations as shown in Fig. [2j 
but separately for each of the gap subsamples, as indicated by 
different symbols. Clearly, groups with different amounts of 
luminosity gap have distinctive Mh(L c ) relations, especially 
at the massive end. For the same central galaxy luminosity, 
groups with a smaller luminosity gap (or larger L 2 ) tends to 
be more massive, especially for groups whose central galaxy 
luminosities are larger than IO 10 0 h~ 2 Lq. For fainter cen¬ 
tral galaxies, the Mh(L c ) relations are similar for the three 
luminosity gaps, due to the small scatter in the halo mass - 
central galaxy luminosity relation. We may notice the much 
larger reduction of error bars at bright central luminosity end 
for large gap groups with log L gap > 1.0 in SAM mock sam¬ 
ple. This is also due to the fact that in SAM the luminosities 
of central galaxies in a large number of low mass groups are 


overestimated compared to the other three mock samples. As 
we mentioned before, this overestimation of luminosity may 
come from the incorrect color modeling of these galaxies. 

The strength of the luminosity-gap dependence increases 
monotonically with the central galaxy luminosity. All these 
suggest that the luminosity gap can be used to reduce the scat¬ 
ter in the Mh — L c relations, particularly at the massive end. 

To model halo mass using also luminosity gap, we formally 
write 

log M h (L c , Lgap) = log M h (L c ) + A log M h (L c , L gap ), 

( 11 ) 

where the first term on the right side is the empirical rela¬ 
tion described by Eq. (fToh . while AlogM/ l (L c ,L gap ) rep¬ 
resents the amount of correction to the M^(L C ) relation with 
the help of L gap . Fig. Q] shows the relation between lumi¬ 
nosity gap and A log Mh for given central galaxy luminos¬ 
ity, where A log Mh is defined by the ratio between the true 
halo mass and that predicted by eq. (ITOt . Shown in different 
panels are for halos (groups) in different L c bins, centered at 
logL c = 9.9, 10.2, 10.5, 10.8, with a bin width equal to 0.1, 
respectively. 

The red solid points represent the median A log Mh for 
each luminosity gap subsamples constructed from the CLF1 
sample, with error bars indicating the 68 percentile. In all 
cases A log Mh decreases with increasing luminosity gap, 
consistent with the trend seen in Fig 0 where one sees that 
halos with a smaller luminosity gap are more massive. The 
overall amplitude of A log Mh increases with the central lu¬ 
minosity, consist with the growing sizes of errorbars at the 
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log [L e / h - 2 L 0 ] log [L c / h _ 2 L 0 ] 

Fig. 5.— Similar to the results shown in Fig. [2] but here we compare the original (open circles) log Mh and luminosity gap pre-corrected (solid circles) 
log M' h halo masses. Obviously, the latter has significantly reduced scatters in log Mh, especially at the massive end. The ratio between these two error bars, <5, 
are shown in the small panels with dashed lines. 



log [M h / h->M 0 ] 


log [M h / h - 1 M 0 ] 


Fi g. 6 .— Comparison to the true halo mass for four mock samples. Red points are halo masses obtained using Eq. [H while black points are obtained using 
Eq.EO Red solid lines show the relation when log Ml equals to log Mh- The standard variance of galaxies from the red solid lines, a, are shown in the small 
panels. The red dotted line corresponds to red points of the above panels, while black one corresponds to black points. 
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large L c end shown in Fig. [2] For different L c the asymptotic 
values of A log Mh are also different. 

We use the following functional form to model 

A log M h {L c , Lgap), 

A log M h (L c , L gap ) = rj a exp(r] b log -^gap) + rjc , ( 12 ) 

where parameters r] a , ij b and r) c may all depend on L c . The 
parameter r/ b in general has a negative value which describes 
the decline rate of A log Mh with log L gap . Parameter -q c de¬ 
scribes the asymptotic value at large log L gap , while param¬ 
eter ij a + r] c is the value of A log at log L gap = 0. To 
gain some idea of the functional forms of these three param¬ 
eters, we first fit to the data for each of the log L c bins shown 
in Fig. |4] and obtain the values of r ) a , r/ b and rj c . According 
to their general dependence on log L c , we use the following 
forms to model these three parameters as functions of log L c : 

Va(Lc) = exp (log L c — Pi) 

Vb(L c ) = a 2 (logL c - /3 2 ), (13) 

V c{L c ) = - (logic -fa) 13 

which in total has five free parameters. We fit all the data 
for CLF1 (as shown partly in Fig. 0 in the luminosity 
range 10.9 > log L c > 9.7 with the above model using a 
Monte-Carlo Markov Chain (MCMC) to explore the likeli¬ 
hood function in the multi-dimensional parameter space (see 
Yan, Madgwick & White 2003; van den Bosch et al. 2005b 
for more detail). Then, we chose the parameter set has the 
highest likelihood (minimum x‘ 2 ) to be the model parame¬ 
ters. The best fit values together with their 68 % confidence 
ranges are listed in the first row of Table [2] This set of pa¬ 
rameter values are valid for central galaxies in the luminosity 
range 10.9 > log L c > 9.7. For central galaxies outside this 
range, one can apply the correction factor at log L c = 9.7 
for fainter galaxies and log L c = 10.9 for brighter galaxies. 
At the faint end, the asymptotic correction factor A log at 
log L c = 9.7 is small. At the bright end, the reason to cut at 
logL c = 10.9 is that our data for massive halos are statisti¬ 
cally poor and not included in the fitting and that the fitting 
function for the Mh(L c ) relation (Eq. IlOi l starts to deviate 
from data at logL c > 10.9. 

The best fit A log as a function of log L gap is 
shown in each panel of Fig. [4] as the solid line. 
Obviously, the model describes the overall behavior in 
A log Mh(L C: i gap ) remarkably good. For comparison, we 
also show A log Mh{L Cl L gap ) as a function of logL gap ob¬ 
tained from the CLF2, SAM and SHAM samples in each 
panel of Fig@] using different symbols as indicated. Al¬ 
though coming from completely different galaxy formation 
models and/or cosmological parameters, the halo mass cor¬ 
rection factors for other samples are quite similar to that for 
the CLF1 sample. Such agreement indicates that our model 
for A log Mh(L c , T gap ) is quite insensitive to the details of 
galaxy formation and cosmology. 

Since the main purpose of this paper is to obtain a better 
halo mass estimator, we compare the amount of scatter in the 
original log Mh(L c ) relation and the new log Mh(L C: L gap ) 
model. To perform such a comparison, we define a ‘pre¬ 
corrected’ halo mass 

log M' h = logM h - Alog M h (L c ,L gap ) , (14) 

and check if the scatter in the M' h (L c ) relation is significantly 
reduced relative to that in the Mh(L c ). If the correction were 
perfect, the scatter in the M' h (L c ) would be reduced to 0. 


Fig. [3 demonstrates the performance of the corrected 
model. Here results for the original Mh(L c ) relations are 
shown as the open circles and the M' h (L c ) relation as the solid 
points. The error-bars are all 68 percentiles around the median 
values in the L c bins. It is clear that the scatter in M’ h (L c ) 
is significantly reduced relative to the original relation, espe¬ 
cially for massive halos/groups where the scatter is reduced 
by about 50%. For clarity, we define the ratio between the 
corrected and original error bars as 5 in the sub-panels. We 
can see, 6 close to 1.0 at low mass ends indicate that two er¬ 
ror bars is about the same. Then S drop to about 0.5 for four 
samples show that the corrected error bar can reduce 50% of 
the original one. Note that the correction model A log Mh is 
calibrated by the CLF1 sample alone, but its application to 
other samples also leads to significant reduction in the scatter. 
These results indicate clearly that adding of even just the sec¬ 
ond brightest galaxy can give a significant improvement of the 
mass estimation over using the central (the brightest) galaxy 
alone. 

In order to further test the reliability of this halo mass esti¬ 
mation, we compared the halo mass log Ml obtained from the 
galaxy luminosity using Eq. QT|to the true halo mass log A//, 
for each group which are directly available in our four mock 
samples: CLF1 CLF2, SAM and SHAM galaxy catalogs. Fig. 
[ 6 ] shows this comparison. The red points represent the halo 
masses estimated using only L c , i.e., using Eq. [TO] While the 
black points represent the halo masses estimated using both 
L c and L gap (Eq. ITO . In order to quantify the scatter with 
respect to the line of equality (log Ml = log Mh) in each 
panel, we measure the standard variance er between estimated 
and true halo mass for each galaxy group. The er is defined 
as: 

ElLi(log M h - log MlY 
n — 1 

where n is the number of groups in each halo mass bin with 
width Alog Mh = 0.1. The results, shown in the small pan¬ 
els, with red dotted lines represent the standard variance for 
only using L c and black ones represent the estimation by us¬ 
ing both L c and L gap , indicate that the estimations for halo 
mass are all improved for four mock samples. Overall, adding 
the second brightest member galaxy can roughly reduce the 
scatter by ~ 0.3dex at the massive end. 

We also noticed that the halo masses estimated using L c 
alone has somewhat systematical deviation from the true halo 
masses at the massive end. Here the systematic deviation is 
mainly induced by the so called Eddington bias where there 
are significantly more low mass halos than massive halos that 
can scattered to larger ones. As we can see, with the help 
of the second brightest member galaxy, this deviation de¬ 
creases substantially. And we will see later that using ad¬ 
ditional member galaxies, such systematic deviation indeed 
disappears finally. 

3.4. Using subsequent member galaxies 

Having demonstrated that the luminosity gap between cen¬ 
tral and the brightest satellite galaxy can be used to better 
constrain the halo mass estimation, we come to subsequent 
member galaxies. Since the richness is also an indicator for 
halo mass and correlated with L c (Shen et al. 2014), we may 
expect the luminosity gap between central and the other satel¬ 
lite galaxies also contains information of halo mass. In this 
subsection, we will extend our study to fainter member satel- 
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TABLE 2 

Parameters for the A log M h function 
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Fig. 7.— Similar to the results shown in Fig. [3 but here for CLF1 mock sample with three, five, seven and nine brightest member galaxies been observed and 
used in calculating the luminosity gap. Groups within three different luminosity gap ranges are shown using different symbols as indicated in the upper and lower 
right panels. Here again, symbols and error bars represent the median and 68 % confidence levels of log M^ in different log L c bins. 


life galaxies. Like before, we still use the ‘luminosity gap’ 
Lgap = L c /L s , but here L s is the luminosity of /-th mem¬ 
ber galaxy in consideration with L s = L3 or L4, etc., where 
L 3 and L,\ are the luminosities of the third and fourth bright¬ 
est galaxies in this group. Note that, in general one can also 
define L s as the total luminosity of all member satellite galax¬ 
ies in consideration. However, as we have tested, this defini¬ 
tion dose not improve any of our results. While the gap thus 
defined has negative values. For simplicity, we adopt L s as 
luminosity of the i-th brightest member galaxy. 


We start with the third brightest member galaxy and stop 
at the ninth. For the fainter member galaxy, the luminosity 
gap will be larger. Overall, the values of L gap are roughly 
in the range 0.0 — 3.0. As before, we split groups into three 
categories according to the value of the luminosity gap. The 
separation criteria for groups with member less than five are: 
small gap (0.5 > logL gap > 0.0), median gap (1.0 > 
log Lg a p > 0.5) and large gap (logL gap > 1.0). The criteria 
for groups with six and more member galaxies are: small gap 
( 1.0 > logL gap > 0 . 0 ), median gap ( 2.0 > logL gap > 1 . 0 ) 
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Fig. 8 .— Similar to the results shown in Fig.[5]but here for four mock samples (CLF1, CLF2, SAM and SHAM from top to bottom rows) with three, five, seven 
and nine brightest members galaxies been observed and used in calculating the luminosity gap. Here we compare the luminosity gap pre-corrected (open circles) 
halo mass log M' h (the second brightest galaxy) in Fig. [5] with the ones obtained from the third, fifth, seventh and ninth brightest member galaxies, respectively. 
Obviously, the latter make further reduction on scatters in log Mi !. especially at the massive end. 


and large gap (log L gap > 2.0). As an illustration, we show 
in Fig. [7] the gap-dependent halo mass distributions, here for 
CLF1 sample only. From top left to right bottom are the re¬ 
sults of luminosity gap between central and the third, fifth, 
seventh and ninth brightest members, respectively. Within 
each panel only halos which contain at least the indicated 
number of member galaxies with log L > 8.0 are used. In 
all four panels, the relationships between central galaxy lumi¬ 
nosity and the halo mass are similar, and the overall luminos¬ 
ity gap dependence is similar to that shown in Fig. [3] for the 
Mh(L c ) relations. Compare to Fig. [3] the separation of halo 
mass between different luminosity gaps increases along with 


the richness. For groups using more than five member galax¬ 
ies, the variance between different gaps are more distinct than 
that by using three or four member galaxies, which suggest 
that a better constraint on the halo mass can be obtained by 
using fainter member galaxies, as we will see below. 

For a given number of members starting from three, we 
first repeat calculations similar to those shown in Fig. Q] and 
then split the sample into several bins according to the central 
galaxy luminosity, starting from log L c = 9.8 to 11.0 with bin 
width equal to 0.1. The A log Mh ( L c , L gap ) — log L gap rela¬ 
tions in the luminosity range 11.0 > log L c > 9.8 are used 
to obtain the best fit parameters that describe the corrections 
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Fig. 9.— Similar to the results shown in Fig. [6] but here we only show the <r values for four samples: CLF1, CLF2, SAM and SHAM. Lines with different 
colors represent different member galaxies that are applied to calculate the luminosity gap. 


defined in Eq. (ITTI) . using the same MCMC algorithm as to the 
central-second cases. The parameters so obtained are listed in 
the second to eighth row of Table|2]for cases using the third, 
fourth, • • •, ninth brightest members, respectively. This set 
of parameters are again only valid for central galaxies within 
the luminosity range used in fitting. For central galaxies with 
luminosities outside this region, one can apply the correction 
factors at log L c = 9.8 for faint galaxies and at log L c = 11.0 
for bright galaxies. We find that the functional form (Eq. 12) 
describes the overall A log Mh(L c , L gap ) remarkably well for 
every central galaxy luminosity bins considered, and for all 
cases regardless which member galaxy used. 

Note that the parameters listed in Tableware obtained from 
the CLF1 sample alone. To see if the same correction model 
can also be applied to other samples, we have checked the 
A log Mh distributions for the CLF2, SAM and SHAM sam¬ 
ples as well, and found that the model also works well. 

To show the performance of the correction model, we use 
the same definition of ‘pre-corrected’ halo mass log M' h as 
described in Eq. (Tufli . To distinguish with the one obtained 
from the second brightest member galaxy, we denote it with 
log M”. As demonstrated in Fig. 5, the scatter is significantly 
reduced in the M’ h (L c ) relation when the second brightest 
galaxies are used in the halo mass estimate. By using subse¬ 
quent fainter members, we expect the model to be improved 
further. In Fig. [8]we compare the ‘pre-corrected’ halo mass as 
a function of L c between cases using different member galax¬ 
ies obtained from the CFF1, CFF2, SAM and SHAM samples 
in different row panels as indicated. The blue open circles 
with error bars show the median and 68 percentile of log M' h , 
while the solid squares with error bars show the same quanti¬ 
ties of log M[[. 


We see that, in all four samples, involving fainter member 
galaxies do make further reduction in the scatter of Mh(L c ) 
relation compared to using the second brightest galaxies. The 
correction method works more effectively for groups with 
larger central galaxy luminosity flog L c > 10.5). Overall 
there is about 40% additional reduction in the scatter at the 
massive end for a fainter (e.g. the fifth - seventh brightest) 
member galaxy. For lower luminosities (10.5 > log L c > 
10 .0), the scatter itself shows that including the second or 
more fainter member galaxies (log M' h v.s. log M£) does not 
change the scatter significantly, because the halo mass is dom¬ 
inated by the central galaxy itself, and the scatter is already 
quite small even without the correction. 

To better assess the quality of our halo mass estimation with 
additional member galaxies, we check the difference between 
the extracted and true halo masses of the groups similar to that 
shown in the small panels of Fig. [6] We show in Fig. [9] the 
standard variance, cr (Eq. ED, of estimated halo mass from 
true halo mass of each group. In each panel lines with differ¬ 
ent colors represent the estimations by using the second, third, 
fourth • • -, ninth brightest members of CFF1, CFF2, SAM and 
SHAM mock samples, respectively. In general, as expected, 
using fainter member galaxies does improve the performance 
of the estimation, especially for massive halos. Interestingly, 
we also noticed that this improvement is not a monotonically 
growth with the number of member galaxies we concerned. 
When the number of member galaxies involved is larger than 
five, the amounts of improvement actually are quite limited 
and become negligible when using more than seven members. 
This feature is caused by the decomposition of our halo mass 
estimation method: the median L c — M h relation together 
with a gap dependent correction. This decomposition has the 
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advantage of shrinking the halo mass estimation error only in 
very poor groups. As we will show in the next section, for 
groups with more than five members, ‘RANK’ method still 
performs better given good survey condition. Hence, we treat 
the seventh member galaxies as the most appropriate option 
in this paper. In future applications, we will use halo mass 
estimation by involving seven member galaxies at most. 

4. COMPARING WITH THE RANKING METHOD 

Having demonstrated the ability of improving the halo mass 
estimations using luminosity gap between the central and 
satellite galaxies, we would like to see if this correction model 
already achieves the similar good performance of estimating 
halo mass using the ranking of Lq (see Fig. Q}. Hereafter we 
refer these two methods as ‘GAP’ and ‘RANK’, respectively. 

4.1. The effect of richness 

The results from last section show the performance of 
‘GAP’ method somehow depends on how many member 
galaxies are involved. Meanwhile, this difference of richness 
certainly will impact on the performance of ‘RANK’ method 
which is based on the characteristic group luminosity of mem¬ 
ber galaxies. Therefore, it would be interesting to compare 
the performances of ‘GAP’ and ‘RANK’ methods on group 
systems of same richnesses. To make fair comparisons be¬ 
tween groups of the same richness, we update the characteris¬ 
tic group luminosity Lq with L G , which is defined as the total 
luminosity of the brightest member galaxies to be included. 
For example, when we only involve three member galaxies, 
then L' g = Li + L 2 + L 3 . Then, the estimation of halo mass 
is based on the ranking of L' G instead of Lq- 

Fi all 01 shows the results of such a comparison. Similar to 
Fig. [ 6 ]and Fig. [G] or is the standard variance of estimated halo 
masses from true halo masses. In each panel, the solid black 
line is given by the ‘GAP’ method and is the same as that 
shown in Fig. [9] which uses the second, third, fifth and sev¬ 
enth brightest members for CFF1, CFF2, SAM and SHAM 
mock samples, respectively. The dotted red line in each panel 
is given by the ‘RANK’ method using only two, three, five 
and seven brightest member galaxies. 

In general, the performance of the ‘GAP’ method is some¬ 
what better when only a few (less than five) brightest member 
galaxies are considered. Then, this advantage is gradually lost 
when involving more members. In the row only five member 
galaxies are considered, the performance of ‘GAP’ and cor¬ 
rected ‘RANK’ method is about the same. When the number 
of member galaxies reach to seven, the ‘RANK’ method tends 
to be slightly better especially for the SAM sample. Mean¬ 
while, we could also see the impact of including different rich¬ 
ness for each galaxy group on ‘RANK’ method is larger than 
‘GAP’ method. Thus, ‘GAP’ method would be a better option 
when only a few galaxy members can be obtained, especially 
for groups with members less than five. 

4.2. The effect of magnitude limit 

In real observations, one can only observe galaxies brighter 
than the survey magnitude limit. This limiting magnitude in¬ 
turn allows us to have a complete observation of galaxies of 
given absolute magnitude to certain redshift. Here we make 
use of CLF1 mock sample to assess the performances of the 
‘RANK’ and ‘GAP’ methods under different absolute mag¬ 
nitude cuts. Note that, for ‘GAP’ method, there are several 
ways to calculate luminosity gap depending on which mem¬ 
ber galaxy are used. As we found in previous section that the 


‘GAP’ method roughly reaches the best performance at 7th 
brightest member galaxy, therefore, in what follows, when we 
refer to ‘GAP’ method, we mean using the maximum avail¬ 
able up to 7th member galaxies to estimate the halo masses. 

Fig. fTTI shows the comparisons between true halo masses 
directly obtained from the mock samples and the ones esti¬ 
mated using ‘GAP’ (left panels) and ‘RANK’ (right panels) 
methods, respectively. We chose the absolute magnitude cuts 
to be 01 M r — 5 log/i = —19.5 (top) and 01 M r — 5 log/i = 
—21.0 (bottom). Using galaxies brighter than these absolute 
magnitudes in a group, we estimate the halo masses of the 
groups according to their characteristic group luminosity and 
luminosity gap respectively. In general, the performances of 
‘RANK’ is slightly better than ‘GAP’ method, especially with 
a fainter absolute magnitude cut. Once we go to brighter ab¬ 
solute magnitude cut, both the ‘RANK’ and ‘GAP’ methods 
have worse accuracy of halo mass estimation. The amount of 
variance between the true halo masses and the estimated ones 
increase by about 0.1 dex from 01 M r — 5 log h = —19.5 to 
01 M r — 5 log h = —21.0 for both estimation methods. Nev¬ 
ertheless, we see that the performance of ‘GAP’ method is 
approaching to that of the ‘RANK’ method. Here, compare to 
results show in Figj 6 ] we do see that the systematic deviation 
in mass estimation disappears. 

To reveal in more detail the impact of different magni¬ 
tude cuts, we plot the comparison of a for both ‘RANK’ and 
‘GAP’ estimation methods in FigJT2] to much brighter abso¬ 
lute magnitude cuts: 01 M r — 5 log h = —19.5, -20.0, -21.0 
and -22.0. These absolute magnitude cuts roughly correspond 
to redshift completeness limit z = 0.09, 0.103, 0.157, 0.22 in 
the SDSS observation with apparent magnitude limit r = 17.6 
(e.g. Wang et al. 2007). The black lines in Fig{T2] repre¬ 
sent the results given by ‘GAP’ method, while red lines are 
given by ’RANK’ method. Apparently, missing faint mem¬ 
ber galaxies influences both methods. Compare to ‘GAP’ 
method, ‘RANK’ method is more sensitive to the increasing 
of luminosity threshold. In 01 M r — 51og/i = —19.5 panel, 
‘RANK’ method has an obvious advantage over a wide halo 
mass range, while the situation is inverted at 01 M r —5 log h = 
—22.0. These results may indicate that, for very high redshift 
or very shallow observations where the luminosity threshold 
is high, ‘GAP’ method may be a better option. In addition, 
the ‘GAP’ method is not suffered from the volume calcula¬ 
tion which is required in the ‘RANK’ method, thus can be 
easily implemented to surveys with poor geometry and are 
flux limited. 

5. SUMMARY 

Galaxies are thought to form and reside within code dark 
matter halos. Different galaxy formation processes, e.g., star 
forming, quenching, AGN feedback, etc., may occur or dom¬ 
inate in halos of different masses. On the other hand, the ha¬ 
los are the building block of the cosmic web, one can use 
the halo mass functions and clustering properties to under¬ 
stand the nature of our Universe. Within these studies, halo 
mass estimation from various observations is one of the key 
challenges. However, such a mission is difficult, especially 
for poor galaxy systems where only a few brightest member 
galaxies are observed. Within the galaxy formation frame¬ 
work, the central (brightest) galaxies in the dark matter halos 
have a rough scaling relation with the masses of the halos: 
brighter galaxies live in more massive halos. However, the 
scatter within that scaling relation is very large, especially in 
massive halos. In this study, we make use of four mock galaxy 
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Fig. 10.— The comparison of <j for two halo mass estimation methods which involve two, three, five and seven brightest member galaxies in each galaxy group 
of CLF1, CLF2, SAM and SHAM mock samples, respectively. The black solid lines repr esen t the ‘GAP’ method, which are the same as those shown in Fig. |9] 
Red dotted lines represent the results for the ‘RANK’ method which are described in the ED 


catalogues, one based on the SAM, one based on the SHAM 
and the other two on the CLF, to investigate the central-halo 
scaling relation and to find out a way that can significantly re¬ 
duce its scatter. Based on these four mock samples, we probed 
the impact of the luminosity gap, which is defined as the lumi¬ 
nosity ratio between the brightest and other member galaxies 
in the same dark matter halo. In this paper, we take into ac¬ 
count a maximum of nine member galaxies in the modeling. 
We find that the scatter in the central-host halo mass relation 
is luminosity gap dependent, which in turn can be used to re¬ 
duce this intrinsic scatter. The main findings of this paper are 
as follows. 


1. We find that the median halo mass for a given central 
galaxy luminosity can be described by simple relation 
described by Eq. EH however with quite large scatter 
around this median. 

2. The scatter in the halo mass depends both on the central 
galaxy luminosity and the luminosity gap between the 
central and the subsequent brightest member galaxies. 

3. We have obtained a mass correction factor Eq. [12] 
which is independent to the detailed galaxy formation 
models, and thus can be applied to any median halo 
mass - central galaxy luminosity {Mh — L c ) relation to 
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Fig. 1 1.— Similar to Fig. [6] but here the comparisons are carried out between ‘GAP’ and ‘RANK’ methods based on CLF1 mock sample with two different 
absolute magnitude cuts: 01 M r — 5 log h = —19.5 (top) and -21.0 (bottom). 
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Fig. 12.— Similar to Fig. ED but here only show the a values for both ‘GAP’ 
01 M r — 5 log h = —19.5, -20.0, -21.0 and -22.0, respectively. 

get better estimation of the halo masses. 

4. The correction factors can reduce the scatters in halo 
mass estimations in massive halos by about 50% to 70% 
depend on which member (second or seventh) galaxies 
are used. 

5. Comparing this ‘GAP’ method with traditional 
‘RANK’ method, we find that the former performs 
better for groups with less than five members, or 


(red line) and ‘RANK’ (black line) methods for four absolute magnitude cuts: 


in observations with very bright magnitude cut. In 
addition, the ‘GAP’ method does not need to calculate 
the volume in estimating halo masses, and thus is much 
easier to be applied to observations with very small 
volume or with poor geometry. 

The above modeling is very useful for our probing of those 
poor galaxy systems or shallow observations where only a few 
brightest member galaxies in a halo can be observed. Based 
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Fig. 13.— Similar to the results shown in Fig. [9] but here the u values are given by another definition of luminosity gap for CLF1.CLF2. SAM and SHAM 
sample respectively. Lines with different colors represent different number of member galaxies that are applied to calculate the luminosity gap. 


on these limited member galaxies, we can have a fairly good 
estimation of the halo mass, which is important for various 
astrophysical and cosmological studies. 
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APPENDIX 

A DIFFERENT DEFINITION OF LUMINOSITY GAP 

Throughout the paper, we define the luminosity gap to be the difference between the brightest and the i-th brightest member 
galaxies. In general, we can also define the luminosity gap as the difference between the brightest and all satellite galaxies, 
logigap = logL c — log L s , here log L s is defined as the total satellite galaxy luminosity in consideration. For example, if we 
only consider two brightest member galaxies in a group, log L s = £ 2 - If we choose to include four member galaxies in a group, 
then log L s = L 2 + i .3 + L-i, where L 3 and L 4 are the luminosities of the third and fourth brightest galaxies in this group. 

Then, we applied exactly the same procedure as before using this new luminosity gap definition. First, calculate the 
A log Mh{L c , Lgap) — log /-gap relations as in Fig. |4] in the luminosity range 10.9 > log L c > 9.7. Then, use MCMC al¬ 
gorithm to find the best parameter set which can describe the relation defined by eq. [12] and eq. [13] Note that, since we applied 
total luminosity of member galaxies, the luminosity gap log L gap can be negative. Overall, the values of log L gap are roughly in 
the range —1.0 — 2.0. Using the relation obtained by the new definition of luminosity gap, we estimated the halo mass for each 
group in four mock samples and compared these estimated halo mass to the true halo mass using er as before. Fig. [l3]shows the 
performance of this new definition of luminosity gap. Compare to those shown in Fig. |9] we can see the tendency and value of <7 
are quite similar in the two figures. In general, the definition for i-th member galaxy luminosity is slightly better at massive ends. 
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